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(57) This invention provides a frequency estimating 
circuit in which estimating a frequency offset distributed 
outside the observation range based on the frequency 
of a received pilot signal degenerated into the observa- 
tion range. The frequency estimating circuit comprises a 
time-axis to frequency-axis converter circuit for convert- 
ing an input data stream constructed of a sampled 
quantized received pilot signal into a signal on a fre- 
quency axis to determine a power spectrum, and an 



MAX search circuit for searching for the maximum value 
of the power spectrum to determine the frequency f 1 
corresponding to the maximum value of the power spec- 
trum, and a frequency computing circuit for computing 
candidates of a frequency offset distributed outside the 
observation range based on the frequency f-j and for 
outputting them as frequency information f ESX . 
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Description 

BACKGROUND OF THE INVENTION 

FIELD OF THE INVENTION 

The present invention relates to an automatic fre- 
quency control (AFC) for use in a demodulator device 
that demodulates a received signal in a satellite com- 
munication system or satellite mobile communications 
system. 

DESCRIPTION OF THE RELATED ART 

As related initial AFC circuits, a technique of esti- 
mating a frequency offset Af has been proposed based 
on the characteristic of power spectrum that a fre- 
quency corresponding to the frequency offset Af of a 
carrier gives a higher power level when power spectrum 
is derived from a received non-modulated signal. 

Hg. 43 shows an example of the organization of a 
prior art initial AFC circuit, disclosed in TECHNICAL 
REPORT OF IEICE IT90-99 (1990) entitled present 
and Future Perspective in the Research and Develop- 
ment of Digital Modem J (authored by Takeuchi et al., 
published January 23. 1991). 

There are shown a received signal input 1 quan- 
tized and sampled, a low-pass filter (hereinafter referred 
to as LPF) 8 for removing out noise outside the fre- 
quency band of the received signal, a multiplier 101 that 
removes a modulating component by multiplying the 
received signal that is M-phase PSK modulated, Fast 
Fourier Transform (hereinafter refer ed to as FFT) circuit 
2 for converting the received signal expressed in the 
time domain into the frequency domain and for deter- 
mining the power spectrum of the received signal, an 
MAX search circuit 3 for searching the maximum value 
of the power spectrum and determining the frequency 
having the maximum- value power, a frequency divider 

102 for frequency dividing the frequency determined by 
the MAX search circuit 3, a frequency estimating circuit 

103 composed of the multiplier 101, the FFT circuit 2, 
the MAX search circuit 3 and the frequency divider 102, 
and a frequency correction block 7 for removing Af from 
the received frequency based on the frequency informa- 
tion estimated by the frequency estimating circuit 103. 

Figs. 45A and 45B are diagrams showing on the 
frequency axis the result that is obtained by subjecting 
the sampled received signal to discrete Fourier trans- 
form (hereinafter referred to as DFT). 

The operation of the AFC circuit is now discussed. 
In the discussion that follows, for convenience of expla- 
nation, it is assumed that he modulation scheme is M- 
phase PSK, and that the received signal is a baseband 
signal expressed in a complex form. The received signal 
is a discrete signal that is quantized through A/D con- 
version. 

By allowing the received signal 1 quantized and 
sampled at a sampling period of T s to pass through LPF 



8. out-of-band noise is removed. The signal without the 
out-of-band noise is M-multiplied by the multiplier 101. 
M-multiplication removes the modulating component 
even further and the signal becomes a non-modulated 
s signal. 

The signal without noise and modulating compo- 
nent is converted from a time-domain signal into a fre- 
quency-domain signal in the FET circuit 2. By squaring 
the frequency-domain converted signal, power spec- 

10 trum of the received signal is obtained. The FFT circuit 
2 outputs power spectrum within a frequency range 
determined by the sampling period T s . Since the FFT 
circuit 2 performs FFT to the quantized and discrete sig- 
nal to derive power spectrum of the received signal, 

is power spectrum of the received signal is also a discrete 
signal. 

Since the received signal that is time-domain to fre- 
quency-domain converted is a non-modulated signal as 
described above, power spectrum exhibits a peak at a 

20 frequency correlated with the frequency offset Af as 
shown in Fig. 44. Hence Af can be estimated by deter- 
mining the frequency of spectrum that gives a maximum 
power spectrum. The MAX search circuit 3 searches for 
a maximum value of the power spectrum in the power 

25 spectrum that is time-domain to frequency domain con- 
verted by the FFT circuit 2. and determines the fre- 
quency corresponding to the maximum-value spectrum. 

It should be noted, however, that since the fre- 
quency determined by the FFT 2 is already M-multiplied 

30 by the multiplier 101 , the frequency actually observed is 
M times the frequency offset Af . Hence to determine the 
actual frequency offset Af. the frequency from the FFT 2 
is M-frequency divided by the frequency divider 102, 
and thus the actual frequency offset Af is obtained. Fre- 

35 quency correction of the received signal is performed 
based on the frequency information estimated by the 
frequency correction block 7. 

Referring to Figs. 45A and Fig. 45B, DFT is now 
discussed. Fig. 45A is a diagram showing the sampled 

40 received signal. In the sampling of the received signal, 
the received signal is typically over-sampled at a symbol 
rate T times m in view of demodulation. Fig. 45B is a 
diagram showing the signal into which the FFT circuit 2 
time-domain to frequency-domain converts a data 

45 stream for L symbols that are over-sampled by m times. 
Frequency range and frequency intervals observable 
are determined by the Nyquist theorem as follows. 

Frequency range : - f s /2 - f s /2 (1) 

50 

(f s = 1/T s = nVT) 



Frequency interval : 1 /(L • T) 
55 T: Symbol period 

T s : Sampling period 
m: Over-sample count 
L: Symbol count 
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The over-sample count and the symbol count to be 
processed are conventionally set with a view to the fre- 
quency range and frequency interval the entire system 
requires. 

Fig. 46 is a diagram illustrating aliasing characteris- . 
tic of FFT when the received signal is sampled at a sam- 
pling period of T s . Aliasing periodically takes place 
every fs, and when the received signal is subjected to 
DFT, out-of-range components are added (aliasing 
noise is added) within the observable frequency range 
of -fs/2 ~ fs/2. To remove signals outside the observable 
frequency range, the received frequency is passed 
through LPF 8 having a transfer function as shown in 
Fig. 47 prior to being applied to the frequency estimat- 
ing circuit 103 as shown in Fig. 43. The transfer function 
of LPF 8 is designed to remove the components outside 
the observable frequency range. 

SUMMARY OF THE INVENTION 

As described above, since the symbol period T is 
constant, from equation (1) the symbol count L must be 
increased to enhance estimating accuracy in frequency 
estimation, and the over-sample count m must be 
increased to widen the frequency range. Hence to get a 
high-accuracy and wide operational range AFC, the 
amount of data increases, and computational load 
increases accordingly. Depending on types of devices in 
use, the AFC suffers a limit in process capability, and 
fails to meet specified requirements, and further, if the 
frequency offset Af is outside the observable range, no 
frequency estimation is possible. 

The present invention has been developed to 
resolve the above problems, and it is an object of the 
present invention to provide a frequency estimating cir- 
cuit that features a high accuracy level and a wide oper- 
ational range with a small amount of computation 
involved, with a view to aliasing characteristic of DFT. tt 
is another object to provide a frequency estimating cir- 
cuit that features a high accuracy level and a wide oper- 
ational range with a smaller amount of computation 
involved, by decimating the sample value signal in the 
process of DFT 

It is yet another object of the present invention to 
provide an AFC circuit that features a high accuracy 
level and a wide operational range using these estimat- 
ing circuits. 

For achieving the aforementioned object, the 
present invention provides a frequency estimating cir- 
cuit comprising: a frequency detector circuit comprising 
a time-axis to frequency-axis converter circuit for con- 
verting an input data stream that is obtained by sam- 
pling a received signal on a time axis into a power 
spectrum on a frequency axis, and a MAX search circuit 
for searching, in the observation frequency range deter- 
mined by a sampling frequency, for the frequency corre- 
sponding to the maximum value of the power spectrum, 
and a frequency computing circuit for estimating the fre- 
quency offset of the carrier of the received signal 



through computation based on the frequency deter- 
mined by the frequency detector circuit and for output- 
ting frequency information. 

Further, the frequency detector circuit comprises a 

5 plurality of the time-axis to frequency-axis converter cir- 
cuits for converting a plurality of input data streams that 
are obtained by sampling the received signal at a plural- 
ity of different sampling rates into power spectrums on a 
frequency axis in accordance with the respective input 

10 data streams, and a plurality of the MAX search circuits 
corresponding to the plurality of the time-axis to fre- 
quency-axis converter circuits, for determining the fre- 
quencies corresponding to the maximum values of the 
power spectrums converted by the respective time-axis 

is to frequency-axis converter circuits, whereby the fre- 
quency computing circuit estimates the frequency offset 
of the carrier of the received signal through computation 
based on the frequencies determined by the frequency 
detector circuit and for outputting frequency information. 

20 Further, the present invention provides a frequency 
estimating circuit further comprising a plurality of filter 
means having transfer functions that allow a signal 
within the observation frequency range to pass there- 
through, whereby the plurality of the time-axis to fre- 

25 quency-axis converter circuits perform time-axis to 
frequency-axis conversion to the data streams that have 
passed through the filter means. 

Further, the present invention provides a frequency 
estimating circuit further comprising a time difference 

30 circuit for introducing a time difference between the plu- 
rality of input data streams that are obtained by sam- 
pling the received signal at the plurality of different 
sampling rates, whereby the time-axis to frequency-axis 
converter circuits perform time-axis to frequency-axis 

35 conversion to the plurality of input data streams after the 
time difference is introduced by the time difference cir- 
cuit. 

Further, the present invention provides a frequency 
estimating circuit further comprising a decimating circuit 

40 for producing a new data stream by decimating the input 
data stream at predetermined intervals of a (an integer) 
samples, whereby the time-axis to frequency-axis con- 
verter circuit converts the output data stream output by 
the decimating circuit into a power spectrum on a fre- 

45 quency axis, a decimating and frequency detector cir- 
cuit is constituted by the decimating circuit, the time- 
axis to frequency-axis converter circuit, and the MAX 
search circuit, and the frequency computing circuit esti- 
mates the frequency offset of the carrier of the received 

so signal through computation based on the frequency 
determined by the MAX search circuit in the decimating 
and frequency detector circuit and outputs frequency 
information. 

Further, the present invention provides a frequency 
55 estimating circuit further comprising a plurality of the 
decimating and frequency detector circuits for perform- 
ing time-axis to frequency-axis conversion to new data 
streams that are produced by decimating the input data 
stream at different intervals in order to determine power 
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spectrums and for determining frequencies correspond- 
ing to the maximum values of the power spectrums, 
whereby the frequency computing circuit estimates the 
frequency offset of the carrier of the received signal 
through computation based on the frequencies deter- 
mined by the decimating and frequency detector circuits 
and outputs frequency information. 

Further, the present invention provides a frequency 
estimating circuit further comprising a second decimat- 
ing and frequency detector circuit comprising a second 
decimating circuit for producing a new data stream by 
decimating the input data stream at predetermined 
intervals of a (an integer) samples, a second time-axis 
to frequency-axis converter circuit for time-axis to fre- 
quency-axis converting the new data stream into a 
power spectrum on a frequency axis, and a second 
MAX search circuit for searching the maximum value of 
the power spectrum determined by the second time- 
axis to frequency-axis converter circuit and determining 
the frequency corresponding to the maximum value, 
and a second frequency computing circuit for estimating 
the frequency offset of the carrier of the received signal 
through computation based on the frequency deter- 
mined by the second decimating and frequency detec- 
tor circuit and for outputting frequency information, 
whereby the time-axis to frequency-axis converter cir- 
cuit in the frequency detector circuit performs time-axis 
to frequency-axis conversion to the input data stream 
and determines the power spectrum of frequency corre- 
sponding to the frequency information from the second 
frequency computing circuit. 

Further, the present invention provides a frequency 
estimating circuit further comprising a third decimating 
circuit for producing a new data spectrum by extracting 
data from an input data stream so that a resulting sam- 
ple count of the new data stream is 1/b (b is a natural 
number) times the sample count of the input data 
stream, whereby the time-axis to frequency-axis con- 
verter circuit in the frequency detector circuit performs 
time-axis to frequency-axis conversion to the input data 
stream produced by the third decimating circuit, into a 
power spectrum on a frequency axis. 

Further, the present invention provides a frequency 
estimating circuit further comprising a fourth decimating 
circuit for producing a new data stream by extracting 
data from an input data stream so that a resulting sam- 
ple count of the new data stream is 1/b (b is a natural 
number) times the sample count of the input data 
stream, and a second frequency detector circuit com- 
prising a third time-axis to frequency-axis converter cir- 
cuit for time-axis to frequency-axis converting the new 
data stream produced by the fourth decimating circuit 
into a power spectrum on a frequency axis, and a third 
MAX search circuit for searching the maximum value of 
the power spectrum determined by the third time-axis to 
frequency-axis converter circuit and for determining the 
frequency corresponding to the maximum value, 
whereby the time-axis to frequency-axis converter cir- 
cuit performs time-axis to frequency-axis conversion to 



the input data stream into a power spectrum on a fre- 
quency axis having the frequency in the vicinity of the 
frequency information determined by the third MAX 
search circuit in the frequency detector circuit. 

5 Further, the present invention provides a frequency 

estimating circuit further comprising a plurality of the 
second decimating and frequency detector circuits 
which produce new different data streams by decimat- 
ing the input data stream at different intervals, time-axis 

10 to frequency-axis convert the new different data 
streams to determine power spectrums, search maxi- 
mum values of the power spectrums and determine the 
frequencies corresponding to the maximum values, 
whereby the second frequency computing circuit esti- 

is mates the frequency offset of the carrier through com- 
putation based on the frequencies determined by the 
plurality of the decimating and frequency detector cir- 
cuits and outputs frequency information. 

Further, the frequency detector circuit further com- 

20 prises a frequency interpolation circuit for interpolating 
the frequency output by the MAX search circuit. 

Furthermore, the present invention provides an 
AFC circuit comprising a frequency detector circuit 
comprising a time-axis to frequency-axis converter cir- 

25 cuit for converting an input data stream that is obtained 
by sampling a received signal on a time axis into a 
power spectrum on a frequency axis, and a MAX search 
circuit for searching, in the observation frequency range 
determined by a sampling frequency, for the frequency 

30 corresponding to the maximum value of the power spec- 
trum, a frequency computing circuit for estimating the 
frequency offset of the carrier of the received signal 
through computation based on the frequency deter- 
mined by the frequency detector circuit and for output- 

35 ting frequency information, and a frequency correction 
block for producing a plurality of demodulated data 
streams by removing the frequency offset f from the 
received signal based on a plurality pieces of frequency 
information output from the frequency computing circuit 

40 and a determining circuit for determining frequency syn- 
chronization using the plurality of demodulated data 
streams output from the frequency correction block 

BRIEF DESCRIPTION OF THE DRAWINGS 

45 

Fig. 1 is a block diagram showing the frequency 
estimating circuit of the embodiment 1 of the present 
invention. 

Figs. 2A and 2B are diagrams showing the opera- 
so tion of the frequency estimation in consideration of 
aliasing. 

Fig. 3 is a block diagram showing the operation of 
the frequency estimating circuit according to the 
embodiment 1 of the present invention. 
55 Fig. 4 is a block diagram of the frequency estimat- 
ing circuit showing the embodiment 2 of the present 
invention. 

Fig. 5 is a diagram showing the operation of the 
embodiment 2 of the present invention. 
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Fig. 6 is a block diagram of the frequency estimat- 
ing circuit showing the embodiment 3 of the present 
invention. 

Fig. 7 is a block diagram of the frequency estimat- 
ing circuit showing the embodiment 4 of the present 
invention. 

Fig. 8 is a diagram showing an example in which a 
time-difference is introduced in decimation. 

Fig. 9 is a flow diagram showing the operation of the . 
frequency estimating circuit in the embodiment 4 of the 
present invention. 

Fig. 10 is a block diagram of the frequency estimat- 
ing circuit showing the embodiment 5 of the present 
invention. 

Fig. 11 is a comparison diagram comparing the 
DFT results derived from the data stream 1 and the data 
stream 2. 

Fig. 12 is a flow diagram showing the operation of 
the frequency estimating circuit in the embocfiment 5 of 
the present invention. 

Fig. 13 is a diagram showing the method of fre- 
quency estimation that is performed within the observa- 
tion range of the data stream 1 using the DFT result of 
the data stream 2. 

Fig. 14 is a block diagram of the frequency estimat- 
ing circuit showing the embodiment 6 of the present 
invention. 

Fig. 15 is a block diagram of the frequency estimat- 
ing circuit showing the embodiment 7 of the present 
invention. 

Fig. 16 is a flow diagram showing the operation of 
the frequency estimating circuit in the embodiment 7 of 
the present invention. 

Fig. 17 is a block diagram of the frequency estimat- 
ing circuit showing the embodiment 8 of the present 
invention. 

Fig. 18 is a diagram showing the operation of the 
embodiment 8 of the present invention. 

Fig. 19 is a block diagram of the frequency estimat- 
ing circuit showing the embodiment 9 of the present 
invention. 

Fig. 20 is a diagram showing the operation of the 
frequency estimating circuit in the embodiment 9 of the 
present invention. 

Fig. 21 is a flow diagram showing the operation of 
the.frequency estimating circuit in the embodiment 9 of 
the present invention. 

Fig. 22 is a block diagram of the frequency estimat- 
ing circuit showing the embodiment 10 of the present 
invention. 

Fig. 23 is a block diagram of the frequency detector 
circuit (2) in the embodiment 11 of the present inven- 
tion. 

Fig. 24 is a comparison diagram comparing the 
DFT results derived from the data stream 1 and the data 
stream 3. 

Fig. 25 is a block diagram of the frequency detector 
circuit (3) in the embodiment 12 of the present inven- 
tion. 



Fig. 26 is a diagram showing the operation of the 
embodiment 12 of the present invention. 

Fig. 27 is a flow diagram showing the operation of 
the frequency detector circuit (3) in the embodiment 12 
s of the present invention. 

Fig. 28 is a block diagram of the frequency estimat- 
ing circuit showing the embodiment 13 of the present 
invention. 

Fig. 29 is a flow diagram showing the operation of 
10 the frequency estimating circuit in the embodiment 1 3 of 
the present invention. 

Fig. 30 is a block diagram of the frequency estimat- 
ing circuit showing the embodiment 14 of the present 
invention. 

is Fig. 31 is a block diagram of the frequency estimat- 
ing circuit showing the embodiment 15 of the present 
invention. 

Fig. 32 is a block diagram of the frequency estimat- 
ing circuit showing the embodiment 16 of the present 
20 invention. 

Fig. 33 is a block diagram of the frequency detector 
circuit (4) in the embodiment 17 of the present inven- 
tion. 

Fig. 34 is a diagram showing the operation of inter- 
25 potation. 

Fig. 35 is a block diagram of the AFC circuit show- 
ing the embodiment 18 of the present invention. 

Fig. 36 is a flow diagram showing the operation of 
the AFC circuit of the present invention. 
30 Fig. 37 is a block diagram of the frequency correc- 
tion block in the AFC circuit of the present invention. 

Fig. 38 is a diagram showing the operation of the 
frequency correction block in the AFC circuit of the 
present invention. 
35 Fig. 39 is a diagram showing the operation of the 
frequency correction block in the AFC circuit of the 
present invention. 

Ftg. 40 is a diagram showing the operation of the 
frequency correction block in the AFC circuit of the 
40 present invention. 

Fig. 41 is a diagram showing the DFT result with a 
not being a factor of the sample count N. 

Fig. 42 is a block diagram showing the feedback 
type AFC circuit of the present invention. 
45 Fig. 43 is a block diagram of the prior art AFC cir- 
cuit. 

Fig. 44 is a diagram showing the power spectrum of 
non-modulated signal. 

Figs. 45A and 45B are diagrams showing the DFT 
so result of the sampled received signal. 

Fig. 46 is a diagram showing aliasing in Fourier 
transform. 

Fig. 47 is a diagram showing one example of trans- 
fer function. 

55 
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DESCRIPTION OF THE PREFERRED EMBODI- 
MENTS 

Embodiment I. 

5 

One embodiment of the present invention is now 
discussed, in the prior art, the received signal is an M- 
phase PSK modulated signal. When the received signal 
is a modulated signal, the multiplier for removing the 
modulating component and the frequency divider for 10 
dividing the estimated frequency are required as shown 
in Fig. 43, but whether or not these components exist 
does not make change in the discussion that follows. 

In the following discussion, for simplicity of explana- 
tion, it is assumed that the received signal is a non-mod- 15 
ulated pilot signal. For a modulated signal, the multiplier 
and the frequency divider are used as shown in Fig. 43. 
In the following embodiments, to convert the received 
sampled signal from the time domain to the frequency 
domain, DFT is used, as an example. Let f EST represent 20 
frequency information that serves as a candidate for the 
estimated frequency offset (hereinafter referred to as 
Af) of the carrier finally output by the frequency estimat- 
ing circuit. Furthermore, the distribution range of Af is 

25 

f L =i Af sf y 

f L : Lower limit of the distribution range of Af 
f u : Upper limit of the distribution range of Af 

30 

Fig. 1 is a block diagram showing the frequency 
estimating circuit that applies when Af is not necessarily 
within an observable range W1 . Fig. 3 is a flow diagram 
of the operation of the frequency estimating circuit 4. 

There are shown a quantized and sampled input 35 
signal 1 (sample value signal), a DFT circuit 2 1 for con- 
verting the received signal expressed in the time 
domain into the frequency domain and for determining 
the power spectrum of the received signal, an MAX 
search circuit 3 for searching the maximum value of the 40 
power spectrum and for determining the frequency giv- 
ing the maximum value of the spectrum, a frequency 
detector circuit 5 composed of the DFT circuit 2* and the 
MAX search circuit 3, a frequency computing circuit (1) 
60 for estimating the frequency that is a candidate for Af 45 
in view of aliasing based on frequency information f 1 
determined by the frequency detector circuit 5, and a 
frequency estimating circuit 4 comprising the frequency 
computing circuit (1) 60 and the frequency detector cir- 
cuit 5, so 

The operation of the frequency estimating circuit is 
now discussed referring to Figs. 2A and 2B and Fig. 3. 
Now, suppose that Af falls not necessarily within the 
observation range. The DFT circuit 2' converts the sam- 
pled received data stream 1 from time-axis signal to a ss 
discrete signal on a frequency axis (step 102), and 
derives a power spectrum by squaring the converted 
discrete signal (step 103). 

The power spectrum output by the DFT circuit 2* is 



10 

fed to the MAX search circuit 3. The MAX search circuit 
3 searches for the maximum value of the input power 
spectrum in power level within the observation range 
W1 (-fs/2 to +fs/2 in Fig. 2A) (step 104), and determines 
and outputs the frequency f -j of the spectrum resulting in 
the maximum value of the power spectrum (step 105). 

The frequency computing circuit 60 carries out the 
following computation to the frequency f determined by 
the MAX search circuit 3, and outputs an estimate value 
f ES r for Af (step 106). 

f EST =f 1 +m • f s (2) 

m : an integer that meets 

f L *(f 1+ m . f s )sfy 

The frequency computing circuit 60 outputs a plu- 
rality of frequency information f ESX based on the equa- 
tion (2). Although f EST estimated as a candidate for the 
frequency offset is a single in nature, the frequency 
computing circuit 60 is unable to determine m in the 
equation (2), and thus estimates and outputs a plurality 
of candidate frequencies (f t - 2f s , f., - f s , f,, f 1 + f s , etc in 
Fig. 2B). 

Embodiment 2. 

Fig. 4 is a block diagram showing an estimating cir- 
cuit for two input systems of data streams in which the 
received signal is sampled at different sampling rates. 
As shown, designated 61 is a frequency computing cir- 
cuit (2) for estimating Af based on two pieces of fre- 
quency information estimated, and DFT circuits 2a, 2b 
and MAX search circuits 3a, 3b are identical to those 
described in connection with the embodiment 1 , desig- 
nated 5 is a frequency detector circuit comprising DFT 
circuits 2a, 2b and MAX search circuits 3a, 3b. 

TTie operation of the frequency estimating circuit is 
now discussed referring to Fig. 4 and Fig. 5. Let f s1 and 
fs2 represent two sampling rates used in producing data 
streams from the received signal, and let DFT 1 and 
DFT 2 represent the respective results that are obtained 
when DFT circuits 2a, 2b perform DFT to two input data 
streams 1a, 1b sampled. 

Fig. 5 shows the resulting DFT, in which Af is shown 
as degenerated f-n within the observable range deter- 
mined by the sampling rate f s1 and f is shown as degen- 
erated f 12 within the observable range determined by 
the sampling rate f S2 . The MAX search circuits 3a, 3b 
detect maximum power spectrums within the observa- 
ble ranges and then outputs their frequencies as f., 1t f 12 . 

Af is estimated from each of the detected frequen- 
cies f-i 1 . f 1 2 with a view to aliasing as follows. 

f EST = f 11 + n 1 * *si (3) 
f EST= f 12 + n 2 * f S2 (4) 

n 1t r^: integers, each satisfying 
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f L=( f 11 +n l * f SV f 12 +n 2 *fs 2 )" f U 

Therefore, the frequency computing circuit (2) 61 s 
determines f EST by solving equation (5) using the fre- 
quency information f n and f 12 output from the MAX 
search circuits 3a. 3b. 

f 11 +n i * f S1 ^ f 12 + n 2 # f S2 ( 5 ) 

Here again in the same way as in the embodiment 

1 , the frequency estimating circuit is unable to deter- 
mine sole n., , n 2 . and thus outputs a plurality of candi- 
dates. 

Embodiment 3. 

Fig. 6 is a block diagram showing a frequency esti- 
mating circuit having two input systems, in which the 
distribution range of Af is approximately known, and 
LPFs are used to cut off components outside the distri- 
bution range of At As shown, designated 8a and 8b are 
LPFs having a transfer function that allows signals 
within a desired observation range to pass there- 
through, and input signals 1a, lb, a frequency detector 
circuit 5 comprising DFT circuits 2a, 2b and MAX 
search circuits 3a, 3b and a frequency computing circuit 
(2) 61 are identical to those described in connection 
with the embodiment 2. 

The operation of the frequency estimating circuit is 
now discussed. In the same way as in the embodiment 

2, Af is sought within the observable range, and the dis- 
tribution range of Af is within the pass bandwidth of the 
LPF, and the condition of n 1( n 2 in equations (3) and (4) 
is as follows. 

nj , n 2 : integers, each satisfying (|f ES Tl m P ass band- 
width of LPF) 

Therefore, the frequency computing circuit (2) 61 
determines n 1( n 2 that satisfy this condition, estimates 
Af, and outputs an estimate value as the frequency 
information f EST . 

Two LPFs 8a, 8b may have different pass band- 
widths. 

Embodiment 4. 

Fig. 7 is a block diagram showing an estimating cir- 
cuit which estimates frequency by introducing a time dif- 
ference between two input systems. In this case, 
discussed herein is an example of a selector 30 as a 
time difference circuit is incorporated in a frequency 
estimating circuit of the embodiment 2. 

Fig. 8 is a diagram illustrating the operation of the 
circuit with a time difference introduced between two 
input systems, and Fig. 9 is a flow diagram showing the 
operation. Referring to Fig. 7, designated 30 is the 



selector included to introduce 'the time difference 
between the two input systems, and the frequency 
detector circuit 5 and the frequency computing circuit 
(2) 61 remain identical to those described in connection 
with the embodiment 2. 

The operation of the frequency estimating circuit is 
discussed referring to Fig. 8 to Fig. 9. Although, in the 
embodiments 2 and 3, the received signal is sampled at 
different sampling rates, both data processed by the 
DFT circuit 2'a and the MAX search circuit 3a and data 
processed by the DFT circuit 2t> and the MAX search 
circuit 3b are originated at the same time. If Af is not 
time variant, and not varying while data is sampled, a 
time difference may be introduced between two input 
data streams. 

The selector 30 selects the DFT circuit 2*a and 
MAX search circuit 3a (step 110), and f^ is derived 
from a data stream 1a (step 111). 

The selector 30 then selects the DFT circuit 2'b and 
MAX search circuit 3b (step 112), and f 12 is derived 
from a data stream 1b (step 113). Fig. 8 shows one 
example, in which a time difference is introduced in dec- 
imation of input data streams to the DFT circuits in this 
embodiment. 

By introducing the time difference between the 
input data streams, correlation between noises of two 
input data streams is smaller than the correlation that is 
derived from the same received signal, the influence of 
noise is reduced, erroneous detection is reduced, and 
estimation of Af is improved. 

Embodiment 5. 

Fig. 1 0 is a block diagram showing a frequency esti- 
mating circuit that produces a new data stream* by 
extracting data at regular intervals from a quantized and 
sampled data stream and thus estimates frequency. 
Designated 41 is a decimating circuit The decimating 
circuit 41 extracts data, at regular intervals of a samples 
from a sample count N in the input data stream 1, and 
produces a new data stream 1\ Herein, a is a factor of 
the sample count N. The rest of the organization of this 
embodiment is identical to that of the preceding embod- 
iments. 

The organization constructed of the decimating cir- 
cuit 41 , DFT circuit 2' and MAX search circuit is referred 
to as a decimating and frequency detector circuit 51 . 

The operation of the frequency estimating circuit is 
discussed referring to Fig. 1 1 and Fig. 12. 

The decimating circuit (1) 41 extracts data at regu- 
lar intervals of a samples from the sample count N in the 
data stream 1 that is sampled at a sampling period T s , 
thereby producing a new data stream 1 * having a sam- 
pling period of aT s (step 115). 

The DFT circuit 2' converts the new data stream 1 ' 
from time-axis signal to frequency-axis signal (step 
1 16), and squares the frequency-axis signal to derive a 
power spectrum (step 103). The MAX search circuit 3 
searches for the power spectrum having a maximum 
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level and determines the frequency f 2 of the maximum 
power spectrum (step 117). 

Fig. 1 1 shows two cases: in one case, the new data 
stream 1 ' is produced by extracting data on a per a sam- 
ples basis from the data stream 1 and is then subjected 
to DFT, and in the other case, the data stream 1 is 
directly subjected to DFT. 

Now consider N samples in the data stream 1 sam- 
pled at a sampling period T s . Data is extracted on a per 
a samples basis (a is a factor of N) from the data stream 
1 to produce the data stream v having a sampling 
period equal to the original sampling period times a, the 
data stream r is time-axis to frequency-axis converted, 
and the result is compared with the result derived from 
the data stream 1 that is directly time-axis to frequency 
axis converted, showing that the observation range for 
the data stream 1 ' is 1/a times the observation range for 
the data stream 1 . The frequency interval from the fre- 
quency-axis converted data stream 1 ' remains identical 
to that from the frequency-axis converted data stream 1 . 

Fig. 13 illustrates a more specific example with a= 
4. Although through the decimation process, the obser- 
vation range of the data stream V is 1/4 the observation 
range of the data stream 1 as shown, Af is estimated 
within a range equivalent to the observation range of the 
data stream 1 by considering aliasing in the same way 
as in the embodiment 1 . 

Specifically, let f 2 represent the frequency detected 
by the MAX search circuit 3, and f EST is expressed as 
follows. 

f EST = f 2+ m * f s /a (5) 
m: an integer satisfying (f L ^ f 2 + m ■ f s /a = f u ) 

A frequency computing circuit 62 is unable to determine 
m in the equation (5), and thus outputs a plurality of can- 
didates. 

Embodiment 6. 

Fig. 14 is a block diagram showing a frequency esti- 
mating circuit having an LPF, in which the distribution 
range of Af is approximately known, and the LPF is used 
to cut off components outside the distribution range of 
Af. As shown, a data stream 1 , a decimating and fre- 
quency detector circuit 51, and a frequency computing 
circuit (3) 62 remain identical to those described in con- 
nection with the embodiment 5, and the LPF 8 is identi- 
cal to that described in connection with the embodiment 
3. 

The operation of the frequency estimating circuit is 
now discussed. The frequency computing circuit (3) 62 
estimates Af using f 2 in the same way as in the embod- 
iment 5. the condition for m in equation (5) is as follows. 

m: an integer satisfying ( |f EST | s pass bandwidth of 
LPF) 



Thus, m satisfying the above condition is deter- 
mined, and f ES j is output as frequency information. 
Here again, as in the embodiment 1, the frequency 
computing circuit (3) 62 is unable to determine m in 
5 equation (5), and thus outputs a plurality of candidates. 

Embodiment 7. 

Fig. 1 5 is a block diagram of a frequency estimating 

10 circuit, in which the input data stream is decimated by 
two systems of decimating circuit for frequency estima- 
tion: In the figure, an input data stream 1 and a decimat- 
ing and frequency detector circuit 51 are identical to 
those described in connection with the preceding 

is embodiment. Designated 61 is a frequency computing 
circuit that performs computation based on the output 
from the decimating and frequency detector circuit 51 . 

The organization shown here in the figure is 
referred to as a multi-input frequency estimating circuit. 

20 The operation of the estimating circuit is discussed 
referring to Fig. 16. The decimating and frequency 
detector circuit 51 extracts, from the data stream 1 , data 
at intervals p and q, both of which are mutually prime, in 
order to produce two new data streams 2 (hereinafter 

25 referred to as data stream 1 a and data stream 1 b) (step 
121). 

The data stream 1a and data stream 1b are sub- 
jected to DFT, and frequency corresponding to Af within 
each of frequency observable ranges is detected from 
30 the maximum value of the spectrum (step 1 22). 

Let f 2 i and f^ represent the detected frequencies, 
and the estimated value of Af is expressed in consider- 
ation of aliasing as follows. 

35 f E ST = f 2i + n i # f sfc ( 6 ) 

f EST = f 22 +n 2 * f S /c ? 

n-, , n 2 : integers, each satisfying (f L n f EST m i u ) 

40 

Af is estimated by determining n1, n2 that satisfy 
equation (7). 

f 21 +^ • f s 4> = f 22 + n 2 • f s/q ( 7 ) 

45 

The frequency computing circuit (2) 61 solves the 
equation (7) using f21, f22, to estimates Af, and outputs 
an estimated value as frequency information f EST (step 
123). Here again, as in the embodiment I, the frequency 
so estimating circuit is unable to determine n v n 2 in the 
equation (7), and. outputs a plurality of candidates. 

Embodiment 8. 

55 Fig. 1 7 is a block diagram showing a multi-input fre- 
quency estimating circuit 28 having an LPF 8, in which 
the distribution range of Af is approximately known, and 
the LPF 8 is used to cut off components outside the dis- 
tribution range of Af . 
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Hie operation of the multi-input frequency estimat- 
ing circuit is now discussed. The frequency computing 
circuit (2) 61 estimates Af using equation (7) in the same 
manner as in the embodiment 7, and the condition for 
n 1t rig is as follows. 

n 1( n 2 : integers, each satisfying (|fEST| m pass 
bandwidth of LPF) 

Fig. 18 shows one specific example. Data is 
extracted from the data stream 1 at intervals of 3 sam- 
ples and 4 samples to produce two new data streams, 
which are then subjected to DFT. Since frequency 
observable ranges are 1/3 and 1/4 the full data (fre- 
quency resolution remains unchanged), observation 
ranges are |f| ^ 20 Hz, |f| s 15 Hz if the observation 
range W1 is |f| ^ 60 Hz with the pass bandwidth of the 
LPF being the observation range W1. Suppose that 
estimated frequencies are 

f 21 =15 Hz 

f 22 = -5Hz 

then, equation (7) is 

15 + n 1 •40 = -5 + n 2 • 30 

By solving equation (7), the frequency computing circuit 
(2) 61 gives 

n 1= 1 
n 2 = 2 

thus, 

f EST = 55Hz 

Embodiment 9. 

Fig. 1 9 is a block diagram showing the embodiment 
according to the eighth invention of the frequency esti- 
mating circuit, Fig. 20 is a diagram illustrating its opera- 
tion, and Fig. 21 is a flow diagram showing the 
operation. 

In Fig. 19, designated 70 is a secondary frequency 
detector circuit (1) composed of a DFT circuit 2' and an 
MAX search circuit 3, and designated 63 is a frequency 
computing circuit (4) that estimates Af within an obser- 
vation range 1 using the frequency information output 
by a decimating and frequency detector circuit and out- 
puts a plurality pieces of frequency information. 

An input data stream 1, a DFT circuit 2*. an MAX 
search circuit 3 and a frequency computing circuit (1) 60 
are identical to those described in connection with the 
embodiment 1, and the decimating and frequency 
detector circuit 51 is identical to that described in con- 
nection with the embodiment 5. 



The operation of the frequency estimating circuit is 
discussed referring to Fig. 19 to Fig. 21 . The input data 
stream 1 is decimated by the decimating and frequency 
detector circuit 51 at predetermined intervals of a sam- 

5 pies to produce a data stream V (step 115), the data 
stream V is time-axis to frequency axis converted, the 
frequency at the maximum level in a spectrum is deter- 
mined and then output as f 2 (step 1 1 8). Based on f 2 , the 
frequency computing circuit 63 performs computation 

io and outputs frequency information as f ES t (step 1 20). 

There remain candidates within a frequency range 
W1 . Since the range outside the frequency range W1 is 
now going to be observed taking advantage of aliasing, 
the number of the candidates of Af is a x n rf a range 

15 equal to n times the observation range W1 is consid- 
ered (n: an integer). The number of candidates is 
reduced in the following technique. If the data stream 1 
is used, Af exists within the observation range 1 and can 
be estimated. 

20 Thus, a plurality of f EST are first determined in the 
same manner as in the embodiment 2. Fig. 20 shows an 
example with a = 4. The DFT circuit 2' performs DFT to 
the data stream 1 again but at four frequencies only cor- 
responding to Af based on f EST output from the-fre- 

25 quency computing circuit (4) 63 to determine the power 
spectrums. The MAX search circuit 3 searches the max- 
imum one in the four power spectrums; and then out- 
puts the frequency corresponding to the maximum 
power spectrum (second from the leftmost one in Fig. 

30 20) . This spectrum has a frequency of f -i . In this way, the 
use of the secondary frequency detector circuit (1) 70 
reduces the number of candidates to 1/a. 

The frequency computing circuit (1) 60 as a suc- 
ceeding stage performs computation of equation (5) 

35 using f 1 , estimates Af over a range inclusive of the range 
outside the observation range W1 , and outputs an esti- 
mated value as f ES T (step 106). Here again, as already 
described, the frequency estimating circuit is unable to 
determine m in equation (5), and thus outputs e plurality 

40 of candidates. 

Embodiment 10. 

Fig. 22 is a block diagram showing a frequency esti- 
45 mating circuit having an LPF, in which the distribution 
range of Af is approximately known, and the LPF is used 
to cut off components outside the distribution range of 
Af. Fig. 22 is a block diagram of the embodiment 6 of the 
frequency estimating circuit. 
so In the figure, the LPF 8 is identical to that described 
with reference to the embodiment 3, and an input data 
stream 1, a decimating and frequency detector circuit 
51, a frequency computing circuit (1) 60, a frequency 
computing circuit (4) 63 and a secondary frequency 
55 detector circuit (1 ) 70 are identical to those described in 
connection with the embodiment 9. 

The operation of the frequency estimating circuit is 
discussed. Since the distribution range of Af is approxi- 
mately known, the distribution range of Af is within the 
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pass bandwidth of the LPF Therefore, the computation 
of the frequency computing circuit (1) 60 is simply per- 
formed within the pass bandwidth of LPF 8. Specifically, 
the frequency computing circuit (1) 60 computes and 
outputs frequency information f EST with the condition for 5 
m as follows. 

m: an integer satisfying ( |f EST | =i Pass bandwidth of 
LPF) 

10 

Embodiment 1 1 . 

Fig. 23 is a block diagram showing the frequency 
detector circuit (2) according to the embodiment 1 1 of 
the present invention. As shown, designated 44 is a 15 
decimating circuit (2) that extracts data from a data 
stream 1 so that the sample count is 1/b times that of 
the data stream 1 to produce a data stream 1 and des- 
ignated 21 is a frequency detector circuit (2) composed 
of the decimating circuit (2) 44 and a frequency detector 20 
circuit 5. The frequency detector circuit 5 is identical to 
that described in connection with the embodiment 1 . 

Referring to Fig. 24, the operation of the frequency 
detector circuit is discussed. The decimating circuit (2) 
44 extracts data from the data stream 1 consecutively 25 
so that the sample count is 1/b times the data stream 1 
and thus produces the data stream 1". The frequency 
detector circuit 5 performs time-axis to frequency-axis 
conversion to the data stream 1", determines the fre- 
quency of the maximum value power spectrum to result 30 
in Af, and output an estimated value as f 3 . 

Rg. 24 is a comparison diagram that compares the 
result of DFT operation carried out to the data stream 1 
with the result of DFT operation carried out to the data 
stream 1" that is derived by decimating the data stream 35 
1 to for it to have the sample count 1/b times the data 
stream 1 . 

In Fig. 24, (a) shows the data stream 1 having N 
samples with a sampling period of T s . (b) shows the 
data stream 1" that is obtained by consecutively deci- 40 
mating the data stream 1 to be 1/b times the data 
stream 1 . (c) and (d) show the results of DFT operations 
that have been carried out to the data stream 1 and the 
data stream 1 These results are attributed to the char- 
acteristic of DFT. 45 

The DFT result of the data stream 1" is 1/T S times 
the DFT result of the data stream 1 in terms of observa- 
tion range, but the frequency interval of the DFT result 
of the data stream 1" is b/NT s , which is b times that of 
the data stream 1 . 50 

Therefore, although the decimation operation 
reduces the computational load, the frequency detector 
circuit 5 presents a lower frequency estimation accu- 
racy in the DFT result of the data stream 1" than in the 
DFT result of the data stream 1, and the frequency 55 
detector circuit (2) 21 in this embodiment is suitable for 
approximate estimation of frequency. 



Embodiment 12. 

Rg. 25 is a block diagram showing the frequency 
detector circuit according to the embodiment 12 of the 
present invention. There are shown a decimating circuit 
(2) 44, a frequency detector circuit 5, and a secondary 
frequency detector circuit (2) 71 . The organization com- 
prising these component is referred to as a frequency 
detector circuit (3) 22. 

Rg. 26 illustrates the operation of the frequency 
detector circuit, and Fig. 27 is a flow diagram of the 
operation. In Fig. 25, designated 71 is the secondary 
frequency detector circuit (2) 71 composed of the DFT 
circuit 2* and MAX search circuit 3, and the decimating 
circuit (2) 44 and the frequency detector circuit 5 are 
identical to those described in connection with the 
embodiment 11. 

Referring to Fig, 25 to Fig. 27, the operation of the 
frequency detector circuit is discussed. The decimating 
circuit (2) 44 decimates the data stream 1 by 1/b times 
to produce a data stream 1 " in the same manner as in 
the embodiment 1 1 (step 127). The frequency detector 
circuit 5 performs DFT to the data stream 1", deter- 
mines the frequency of the maximum power level spec- 
trum, and outputs frequency information f 3 (step 128). 

The secondary frequency estimating circuit (2) 71 
performs DFT to f 3 that is estimated based on the DFT 
result of the data stream 1" and to the data stream 1 in 
the vicinity of f 3 , and outputs the result of frequency esti- 
mation as frequency information f EST (step 129). 

Fig. 26 shows an example of operation with b=4. 
Since as shown, the frequency interval is quadrupled as 
a result of decimation in the DFT result by the frequency 
detector circuit 5, the secondary frequency detector cir- 
cuit (2) 71 performs afresh frequency estimation in the 
data stream 1 in frequencies in the vicinity of f 3 (a total 
of 9 frequencies f 3 - 1 ~ f 3 + 1 in Fig. 26). In the second- 
ary estimation, the frequency interval becomes quar- 
tered, and the estimation offers a higher accuracy 
compared with the estimation of the data stream 1 

Embodiment 13. 

Rg. 28 is a block diagram of the frequency estimat- 
ing circuit in which the preceding embodiment 12 is 
incorporated in the embodiment 1 . There are shown a 
input data stream 1, a frequency detector circuit (3) 22 
and a frequency computing circuit (1) 60. Rg. 29 is a 
flow diagram showing the operation of the frequency 
estimating circuit. 

The operation of the frequency estimating circuit is 
now discussed. The frequency detector circuit (3) 22 
produces a data stream 1 M by decimating the input data 
stream 1 so that the data stream 1 " has a sample count 
1/b times the sample count N in the input data stream 1 
(step 1 27), the data stream 1 " is then subjected to time- 
axis to frequency-axis conversion, the frequency of 
maximum power level spectrum is determined and thus 
f 3 is determined (step 128). 
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DFT is again performed to only the vicinity of f 3 and 
frequency information is output as f ES j (step 129). The 
frequency computing circuit (1) 60 estimates Af in view 
of aliasing in the same manner as in the embodiment 1 . 
and frequency information f ES T < s estimated (step 130). 

Embodiment 14. 

Fig. 30 is a block diagram showing the frequency 
estimating circuit of one embodiment of the present 
invention. Described herein is the embodiment in which 
above described frequency detector circuit described in 
connection with the embodiment 1 1 or the embodiment 
12 is incorporated in the frequency estimating circuit 
described in connection with the embodiment 10. 

As shown, the input data stream 1 , frequency com- 
puting circuit (1) 60, frequency estimating circuit (4) 63 
and secondary frequency detector circuit (2) 71 are 
identical to those described in connection with the 
embodiment 10, the decimating circuit (1) 41 is identical 
to that described in connection with the embodiment 5, 
and the frequency detector circuit 5 and decimating cir- 
cuit (2) 44 are identical to those described in connection 
with the embodiment 11. 

The operation of this embodiment is discussed 
referring to Fig. 30. By allowing the input data stream 1 
to pass through the decimating circuit (1) 41 and the 
decimating circuit (2) 44, a new data stream is pro- 
duced. When the new data stream is subjected to DFT, 
both characteristic 1 and characteristic 2 of DFT already 
described are applied to the DFT result, and when the 
frequency estimated by the frequency detector circuit 5 
is compared with the DFT result of the input data stream 
1 , the computational load involved is smaller, the esti- 
mated f is degenerated into the observation range W2 
and its discrimination is poor. 

For this reason, the secondary frequency detector 
circuit (2) 71 performs afresh frequency estimation 
using the input data stream 1 in the vicinity of the fre- 
quency estimated by the frequency computing circuit (4) 
63, and then estimates frequency information f ES t The 
order of arrangement is not important in that the deci- 
mating circuit (1) 41 may be arranged ahead of the dec- 
imating circuit (2) 44 and vice versa, or both circuits may 
perform decimation concurrently 

Embodiment 15. 

Fig. 31 is a block diagram showing the frequency 
estimating circuit of one embodiment of the present 
invention. Described herein is the embodiment in which 
the secondary frequency detector circuit described in 
connection with the embodiment 12 is incorporated in 
the frequency estimating circuit described in connection 
with the embodiment 7. 

As shown, a decimating and frequency detector cir- 
cuit 51, a frequency computing circuit (2) 63, and a 
mufti-input frequency estimating circuit 2*8 are identical 
to those described in the embodiment 7, the secondary 



frequency estimating circuit (2) 71 is identical to that 
described in the embodiment 12, and the frequency 
computing circuit (1) 4 is identical to that described in 
the embodiment 4. 

s The operation of the frequency estimating circuit is 

now discussed referring to Fig. 31 . Two decimating and 
frequency detector circuits 51 decimate the input data 
stream 1 at two different intervals which are mutually 
prime to produce data streams 2. Each data stream is 

10 subjected to time-axis to frequency-axis conversion, 
and frequencies in the vicinity of the frequencies of 
maximum power level spectrums in the data spectrums 
2 are output and frequencies corresponding to several 
power spectrums from the power spectrum of highest 

is level, as f 2 i and f 22 , are output. The frequency comput- 
ing circuit 61 computes equation (7) using f 21 , f 22 ,and 
outputs a plurality of frequency information f -j. 

In response to the time-axis to frequency-axis con- 
version of the input data stream I, the secondary fre- 

20 quency estimating circuit (2) 71 searches power 
spectrums in the vicinity of the frequency information f 1 
from the frequency computing circuit 61 to find and then 
output the frequency of the maximum value power spec- 
trum. The frequency computing circuit (1) 60 computes 

25 equation (2) based on the frequency information to esti- 
mate frequency outside the observation range 1 , and 
then outputs the estimated frequency information f EST . 

Embodiment 16. 

30 

Fig. 32 shows the embodiment which is constructed 
of the preceding embodiment with an LPF 8 added, 
wherein its frequency offset is approximately known. 
The operation of this embodiment is now discussed 

35 referring to Fig. 32. The operation of the multi-input fre- 
quency detector circuit 28 and secondary frequency 
detector circuit (2) 71 is identical to that of the preceding 
embodiment. In this embodiment, since the distribution 
range of the carrier is approximately known, the fre- 

40 quency estimating circuit results in an acceptable fre- 
quency estimation within the observation range that is 
determined by a sampling period at which the input data 
stream 1 is sampled. 

The frequency computing circuit (1) computes 

45 equation (2) based on the output information from the 
multi- input frequency detector circuit 28 with the condi- 
tion for m as follows and outputs f EST . 

m: an integer satisfying ( |f ES Tl - Pass bandwidth of 
so LPF) 

Embodiment 1 7. 

Fig. 33 is a block diagram of the frequency detector 
55 circuit having frequency interpolation circuit therewithin. 
Figure shows the embodiment in which a frequency 
interpolation circuit is incorporated in the detector circuit 
of the embodiment 1 . 

The operation of this embodiment is now dis- 



11 



BNSDOCID: <EP 0788225A1J_> 



21 



EP0 788 225 A1 



22 



cussed. The DFT circuit 2* and the MAX search circuit 3 
perform frequency estimation using the power spectrum 
of the input data stream 1 . Although the frequency that 
gives a maximum power spectrum is conventionally out- 
put as frequency information ^ , the frequency that can 5 
be observed as a result of DFT is discrete value, and 
estimation error takes place if Af fails to agree with the 
frequency point of observation. 

For this reason, to estimate frequency at a higher 
accuracy level, the output frequency of the MAX search 10 
circuit 3 derived from the maximum power spectrum is 
interpolated. As an example of interpolation, Lagrange's 
secondary interpolation is used herein. Fig. 34 shows 
the operation of interpolation. In Fig. 34, it is assumed 
that f t is at a maximum. 75 

A frequency interpolation circuit 20 interpolates 
between three power values, a power value at f-| and 
power values at two points (fi - 1 , fi + 1) adjacent to f 1t 
and then estimates the frequency that is supposed to 
give the true peak power value. The frequency interpo- 20 
lation circuit 20 performs interpolation in this way, and 
outputs afresh as f ^ the frequency resulting from inter- 
polation. 

Embodiment 18. 25 

Fig. 35 is a block diagram of an AFC circuit as the 
embodiment 1 1 according to the eleventh invention in 
which the frequency estimating circuit is incorporated. 
As shown, an input data stream 1 and a frequency esti- 30 
mating circuit 4 are identical to those described in con- 
nection with the embodiment 1. There are also shown a 
frequency correction block 7 for correcting a received 
signal and a determining circuit 2*9 for determining syn- 
chronization. 35 

The operation of this embodiment is now discussed 
referring to the flow diagram in Fig. 36. 

The frequency estimating circuit 4 time-axis to fre- 
quency-axis converts the input data stream 1 to deter- 
mine a power spectrum, and determines the frequency 40 
corresponding to the maximum power, computes a pre- 
determined equation (equation (2) in this case) to esti- 
mate Af, and outputs a plurality of frequency information 
f est (step 1 32). 

The frequency correction block 7 removes Af from as 
the input data stream 1 using the plurality of frequency 
information f EST (step 133) to produce a plurality of 
demodulated data streams. The determining circuit 2'9 
determines the synchronization of frequency using the 
plurality of demodulated data streams output from the so 
frequency correction block (step 134), and outputs its 
determination result to the frequency correction block. 

The frequency correction block 7 is now discussed 
further in detail. Fig. 37 shows an example of the organ- 
ization of the frequency correction block 7. f ES T - A6 ss 
converter 81 determines Af based on f £ST output by the 
frequency estimating circuit 4. A reversed phase gener- 
ator circuit 82 determines reversed phase angle of rota- 
tion at each sampling point using Af determined by the 



f EST - A6 converter 81 . 

A phase reverse rotating circuit 83 removes a quan- 
tity of phase change from the received data stream 
based on the input data stream 1 and the angle of rota- 
tion corresponding to each data. Af is removed in this 
way. As data to be demodulated, the received signal 
having no Af is output to the determining circuit 29 and 
a demodulator block at the next stage. 

Next, Af is removed from the received signal 1 
based on a plurality of frequency information f EST to 
produce a plurality of data streams. One method of 
removing Af is discussed next. 

Since the received signal is a pilot signal (non-mod- 
ulated signal), the phase of the received signal is 0 with 
Af= 0 regardless of sampling timing (T s ). If Af* 0, how- 
ever, the phase of the received signal varies with time 
as shown in Fig. 39, and the quantity of phase change 
A0 is 

A6 = Af x 360 x T s 

AG : Phase change per sample (degrees) 
Af : Frequency offset estimated (Hz) 
T s : Sampling period (seconds) 

Since the phase change is determined by determin- 
ing using the frequency information f EST estimated by 
the frequency estimating circuit and the sampling period 
T s , the phase at each symbol becomes 0 if the phase 
change is fed back at each symbol. The frequency cor- 
rection block 7 thus determines A6 based on f EST esti- 
mated by the frequency estimating circuit 4. A complex 
multiplier is used to correct the phase at each sampling 
point as shown in Fig. 40. The data stream thus cor- 
rected is output to the later stage. 

One example of determining method is now dis- 
cussed. A known pattern detection such as UW detec- 
tion is performed. Since the distribution range of Af is 
generally finite (f L ^ Af § f y ), m in equation (2) is a finite 
number. Suppose that the number of a plurality ol fre- 
quency information output by the frequency estimating 
circuit 4 is K, and that these pieces of frequency infor- 
mation are f esti - *estk- 

The frequency correction block 7 produces a 
demodulated data stream by removing Af correspond- 
ing to f E sn f rom * ne received data, and then outputs it 
to the determining circuit 29. The determining circuit 29 
detects a known pattern used in the system such as UW 
in the demodulated data stream. When no known pat- 
tern is detected, the determining circuit 29 outputs no- 
detection signal to the frequency correction block 7. 

Receiving the no-detection signal, the frequency 
correction block 7 outputs a demodulated data stream 
that is produced by removing Af corresponding to next 
frequency information f E sT2 from the received signal. 
This process is repeated until a known pattern is 
detected. 

When the known pattern is detected, the determin- 
ing circuit 29 outputs a detected signal to the frequency 
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correction block 7. When the frequency correction block 
7 receives the detected signal, it fixes a frequency to be 
corrected, with a frequency corresponding to it and cor- 
rects with the frequency thereafter. If no known pattern 
is detected in connection of all K candidates, the AFC 
operation is regarded as failed. 

The prior art AFC circuit can estimate Af only within 
the observable range. Thus, the prior art AFC suffers 
the disadvantage that it cannot perform frequency esti- 
mation when Af is large and outside the observable 
range even with no influence of noise. According to the 
present invention, however, the observable range is infi- 
nite, and Af is estimated without fail with no influence of 
noise or the like, and a highly accurate and reliable syn- 
chronization is established. 

Alternative embodiment 1 

In the embodiment 9 and embodiment 10, the data 
stream 1 ' is used to estimate Af within the observation 
range W2, and a plurality of candidates for Af within the 
observation range W1 are estimated by the frequency 
computing circuit (1) 60, DFT is performed only these 
plurality of frequencies by the secondary frequency esti- 
mating circuit (1) 70, and f EST within the observation 
range W1 is determined, and the use of the characteris- 
tic 1 of DFT is on condition that a must be a factor of 
sample count N. 

Described herein is the embodiment in which a is 
not a factor of the sarnple.count N. 

Fig. 41 shows the DFT result when a is not a factor 
of the sample count N. The operation of the embodi- 
ment is now discussed referring to Fig. 41 . When a is 
not a factor of the sample count N, N cannot be divided 
by a, and if data extraction is performed from the first 
sample at intervals of a samples, d data remains at the 
end of the data stream (d < a), an observation time 
length for DFT is shortened, and the observation time 
length is expressed by equation (8). 

N* • T s = (N-d) • T s (8) 

As a result, the frequency interval becomes longer, 
and it can be expressed by equation (9) as follows. 

1/(N' • T s ) = 1/{(N-d) • T s } (9) 

In this way, since the frequency point of observation 
with the data stream V subjected to DFT is different 
from the frequency point of observation with the data 
stream 1 subjected to DFT, the secondary frequency 
estimating circuit (1) 70 cannot observe the frequency 
that corresponds to the estimated value f E sT with the 
data stream V. In such a case, therefore, the secondary 
frequency estimating circuit (1) 70 performs DFT to the 
data stream 1 at a plurality of points of observation in 
the vicinity of f EST . 



Alternative embodiment 2 * 

In the embodiment 7 and embodiment 8, as well, 
the use of the characteristic 1 of DFT is on condition that 
s a is a factor of the sample count N. 

Described herein is the embodiment in which a is 
not a factor of the sample count N. When a is not a fac- 
tor of the sample count N in both decimating circuits 41 
or one of decimating circuits 41 , the frequency interval 
10 is expressed by equation (9), estimated values f 21> I22 
from two decimating and frequency detector circuits are 
ones over different frequency intervals, and thus n 1( 
satisfying equation (5) do not exist. 

In such a case, the frequency computing circuit (2) 
15 63 uses equation (10) modified from equation (5) to 
determine n 1f n 2 . 

f 21 + n 1 • f/p = f 22 + n 2 • f/q (10) 

20 Alternative embodiment 3 

In the embodiment 1 8, the AFC is arranged in the 
form of feedforward type, but alternatively, it may be 
arranged in the form of feedback type. 

25 Such an arrangement is shown in Fig. 42. 

In Fig. 42. a mixer 18 multiplies the output from an 
oscillator for down-conversion (hereinafter referred to as 
D/C) by a received signal input to remove Af from the 
received signal, and a D/C oscillator 19 outputs the fre- 

30 quency corresponding to the frequency information f EST 
output by a frequency estimating circuit 4, and the fre- 
quency estimating circuit 4 and a determining circuit 29 
are identical to those described in connection with the 
embodiment 18. 

35 In the same way as in the embodiment 18, the fre- 
quency estimating circuit 4 estimates Af, and outputs a 
plurality of frequency information. The D/C oscillator 19 
outputs the frequency corresponding to the frequency 
information, and the mixer 18 multiplies the output from 

40 the D/C oscillator 19 by the input received signal to 
remove Af from the received signal. The determining cir- 
cuit 29 determines the synchronization of frequency 
using the demodulated data stream output by the mixer 
18, and outputs its determination result. 

45 Thereafter, the determination of the synchroniza- 
tion of frequency are performed in the same way as in 
the embodiment 18. 

Alternative embodiment 4 

50 

In the frequency estimation in the embodiment 2 - 
embodiment 4, embodiment 7, embodiment 8. embodi- 
ment 1 5 and embodiment 1 6. two processings are per- 
formed, and increasing the number of processings may 
55 reduce detection error and allow Af to be estimated at a 
good accuracy level. 

Although in these embodiments, processings are 
parallel, a single decimating and frequency estimating 
circuit may be allowed to perform a plurality of process- 
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ings using buffers if processing time presents no prob- 
lem, in this way, the decimating and frequency 
estimating circuit may be of a single arrangement, and 
its construction is simplified. 

Alternative embodiment 5 

The LPF is provided to remove components outside 
the band of the signal in the embodiment 3, embodi- 
ment 6, embodiment 8, embodiment 10. embodiment 
14 and embodiment 16, but such an LPF is needless if 
the aliasing noise and noise within the desired observa- 
tion range are negligibly low (with an excellent S/N 
ratio). Such an arrangement simplifies the construction 
of each of the above embodiments. 

Alternative embodiment 6 

In the embodiment 18, the frequency correction 
process and the determining process are repeated to 
f ESTi ~ *estk unti ' a known pattern is detected, and if 
the distribution of frequency offset Af has any particular 
feature, the processes may be carried out in the order 
that is suited to the feature. 

For example, the distribution of Af is in the vicinity of 
0. the processing may be first performed to frequencies 
corresponding to f EST close to 0, and then f EST may be 
selected in the order of farther away from 0 for process- 
ing. In this way, the time to search for Af is shortened, 
and thus the time to establish synchronization is accord- 
ingly shortened. 

Alternative embodiment 7 

In the embodiment 18, once a known pattern is 
detected, the remaining candidates are not subjected to 
the frequency correction process and the determining 
process. Since there is a possibility of detection error 
under the influence of noise, the frequency correction 
process and the determining process may be performed 
to the remaining candidates even if the known pattern is 
detected at a frequency. 

When the known pattern is found in a single candi- 
date only, the frequency correction block 7, upon receiv- 
ing the detected signal, fixes a frequency to be 
corrected, with a frequency corresponding to it, and cor- 
rects thereafter with the frequency. If the known pattern 
is found in a plurality of candidates, the determination of 
the frequency may be left to the demodulator block at 
the succeeding stage. In this way, detection error due to 
noise or the like is prevented, and Af is estimated at an 
improved accuracy. 

Although this embodiment has been described in 
connection with the AFC circuit into which the embodi- 
ment 1 is incorporated, detection error due to noise or 
the like may be controlled if the frequency estimating cir- 
cuit in the embodiments 2 through 17 is employed, and 
Af is equally estimated at an improved accuracy. 

As described above, according to the present 



invention, time-axis to frequency-axis conversion is per- 
formed and frequency information is estimated based 
on a plurality of candidates of Af from the frequencies of 
the maximum power spectrum, and thus, the frequen- 

5 cies in a conventionally unobservable range (outside 
the observable range) can be observed, and a wider fre- 
quency range is observed with a data quantity remain- 
ing unchanged from that in the prior art system. 

Frequency estimation is performed in different 

10 observation ranges by using a plurality of data streams 
sampled at different sampling rates, and thus Af in a 
conventionally observable range is estimated within the 
observation range with a smaller quantity of data by tak- 
ing into consideration aliasing in a plurality of frequency 

15 information, and furthermore, with aliasing considered, 
a wider frequency range is observed with the data quan- 
tity remaining unchanged from that in the prior system. 

Since, in the course of frequency estimation, a time 
difference is introduced between two data streams to be 

20 used in frequency estimation, correlation of noise 
between the two data streams is small compared with 
that resulting from the decimation of the same data 
stream I, and the influence of noise is reduced, the 
occurrence of detection error is lessened, and a more 

25 accurate Af is thus estimated. 

By allowing the decimating circuit to decimate the 
input data stream in the frequency estimation, Af is esti- 
mated within the observation range with a smaller quan- 
tity of data than that in the prior art system, and further 

30 the same advantage as the first invention is obtained- 
taking into consideration aliasing. 

As described above, a plurality of decimating cir- 
cuits are used to performed frequency estimation in dif- 
ferent frequency ranges and aliasing is considered in a 

35 plurality of frequency information, and thus f is esti- 
mated within the observation range with a smaller quan- 
tity of data than that in the prior art system, and further 
the same advantage as the first invention is obtained 
taking into consideration aliasing. 

40 Fathermore, as described above, the same advan- 
tage as the fourth invention is provided; and, further all 
data is used to perform afresh frequency estimation and 
thus a single frequency candidate is obtained. In view of 
aliasing, the same advantage as the first invention is 

45 provided. 

Moreover, as described above, the use of the deci- 
mating circuit allows rough frequency estimation with a 
smaller data quantity. 

Frequency estimation is performed with a smaller 

so data quantity after the data stream is decimated, all data 
is afresh subjected to frequency estimation, and thus a 
high discrimination frequency estimation is performed 
wfth a small processing quantity. 

Since the frequency detection result from the fre- 

55 quency estimating circuit is frequency interpolated, fre- 
quency estimation is performed at a high accuracy level. 

In the AFC circuit that demodulates the received 
signal, the input data stream, a plurality of estimatere 
used to determine synchronization, and frequency esti- 
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mation is performed in a wide range and a highly accu- 
rate and reliable synchronization is established. 

Claims 

1 . A frequency estimating circuit comprising: 

a frequency detector circuit (5) comprising a 
time-axis to frequency-axis converter circuit (2) 
for converting a data stream (1 ) that is obtained 
by sampling a received signal on a time axis 
into a power spectrum on a frequency axis, and 
a MAX search circuit (3) for searching, in the 
observation frequency range determined by a 
sampling frequency, for the frequency corre- 
sponding to the maximum value of the power 
spectrum, and 

a frequency computing circuit (60) for estimat- 
ing the frequency offset of the carrier of the 
received signal through computation based on 
the frequency determined by the frequency 
detector circuit and for outputting frequency 
information. 

2. The frequency estimating circuit according to Claim 
1 comprising: 



introducing a time difference between the plurality 
of input data streams that are obtained by sampling 
the received signaJ at the plurality of different sam- 
pling rates, whereby the time-axis to frequency-axis 
5 converter circuits (2a, 2b) per form time-axis to fre- 

quency-axis conversion to the plurality of input data 
streams after the time difference is introduced by 
the time difference circuit. 

w 5. The frequency estimating circuit according to Claim 
1 further comprising a decimating circuit (41) for 
producing a new data stream by decimating the 
input data stream at predetermined intervals of a 
(an integer) samples, whereby the time-axis to fre- 

15 quency-axis converter circuit (2) converts the out- 
put data stream output by the decimating circuit into 
a power spectrum on a frequency axis, a decimat- 
ing and frequency detector circuit (51 ) is constituted 
by the decimating circuit (41), the time-axis to fre- 

20 quency-axis converter circuit (2), and the MAX 
search circuit (3), and said frequency computing 
circuit (62) estimates the frequency offset of the 
carrier of the received signal through computation 
based on the frequency determined by the MAX 

25 search circuit in the decimating and frequency 
detector circuit and outputs frequency information. 



said frequency detector circuit (5) comprising a 
plurality of the time-axis to frequency-axis con- 
verter circuits (2a,2b) for converting a plurality 30 
of input data streams that are obtained by sam- 
pling the received signal at a plurality of differ- 
ent sampling rates into power spectrums on a 
frequency axis in accordance with the respec- 
tive input data streams (1a, 1b), and a plurality 35 
of the MAX search circuits (3a,3b) correspond- 
ing to the plurality of the time-axis to frequency- 
axis converter circuits, for determining the fre- 
quencies corresponding to the maximum val- 
ues of the power spectrums converted by the 40 
respective time-axis to frequency-axis con- 
verter circuits, whereby said frequency com- 
puting circuit estimates the frequency offset of 
the carrier of the received signal through com- 
putation based on the frequencies determined 45 
by the frequency detector circuit and for output- 
ting frequency information. 

3. The frequency estimating circuit according to Claim 

2 further comprising a plurality of filter, means so 
(8a, 8b) having transfer functions that allow a signal 
within the observation frequency range to pass 
therethrough, whereby the plurality of the time-axis 
to frequency-axis converter circuits (2a.2b) perform 
time-axis to frequency-axis conversion to the data ss 
streams that have passed through the fitter means. 

4. The frequency estimating circuit according to Claim 
2 further comprising a time difference circuit (30) for 



6. The frequency estimating circuit according to Claim 
1 further comprising: 

a second decimating and frequency detector 
circuit (51) comprising a second decimating cir- 
cuit (41) for producing a new data stream by 
decimating the input data stream at predeter- 
mined intervals of a (an integer) samples, a 
second time-axis to frequency-axis converter 
circuit (2) for time-axis to frequency-axis con- 
verting the new data stream into a power spec- 
trum on a frequency axis, and a second MAX 
search circuit (3) for searching the maximum 
value of the power spectrum determined by the 
second time-axis to frequency-axis converter 
circuit and determining the frequency corre- 
sponding to the maximum value, and 
a second frequency computing circuit (63) for 
estimating the frequency offset of the carrier of 
the received signal through computation based 
on the frequency determined by the second 
decimating and frequency detector circuit and 
for outputting frequency information, whereby 
the time-axis to frequency-axis converter circuit 
(2) in the frequency detector circuit (70) per- 
forms time-axis to frequency-axis conversion to 
the input data stream and determines the 
power spectrum of frequency corresponding to 
the frequency information from the second fre- 
quency computing circuit (63). 

7. The frequency estimating circuit according to Claim 
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9. 



1 further comprising a third decimating circuit (44) 
for producing a new data spectrum by extracting 
data from an input data stream so that a resulting 
sample count of the new data stream is 1/b (b is a 
natural number) times the sample count of the input 
data stream, whereby the time-axis to frequency- 
axis converter circuit (2) in the frequency detector 
circuit (5) performs time-axis to frequency-axis con- 
version to the input data stream produced by the 
third decimating circuit (44), into a power spectrum 
on a frequency axis. 

The frequency estimating circuit according to Claim 
1 further comprising: 

a fourth decimating circuit (44) for producing a 
new data stream by extracting data from an 
input data stream so that a resulting sample 
count of the new data stream is 1 A) (b is a nat- 
ural number) times the sample count of the 
input data stream, and 

a second frequency detector circuit (5) com- 
prising a third time-axis to frequency-axis con- 
verter circuit (2) for time-axis to frequency-axis 
converting the new data stream produced by 
the fourth decimating circuit into a power spec- 
trum on a frequency axis, and a third MAX 
search circuit (3) for searching the maximum 
value of the power spectrum determined by the 
third time-axis to frequency-axis converter cir- 
cuit and for determining the frequency corre- 
sponding to the maximum value, whereby the 
time-axis to frequency-axis converter circuit (2) 
performs time-axis to frequency-axis conver- 
sion to the ipput data stream into a power spec- 
trum on a frequency axis having the frequency 
in the vicinity of the frequency information 
determined by the third MAX search circuit in 
the frequency detector circuit. 

The frequency estimating circuit according to Claim 
1 , wherein the frequency detector circuit (5) further 
comprises a frequency interpolation circuit (20) for 
interpolating the frequency output by the MAX 
search circuit (3). 
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ing the frequency offset of the carrier of the 
received signal through computation based on 
the frequency determined by the frequency 
detector circuit and for outputting frequency 
information, 

a frequency correction block (7) for producing a 
plurality of demodulated data streams by 
removing the frequency offset f from the 
received signal based on a plurality pieces of 
frequency information output from the fre- 
quency computing circuit, and 
a determining circuit (29) for determining fre- 
quency synchronization using the plurality of 
demodulated data streams output from the fre- 
quency correction block. 



10. An AFC circuit comprising: 

a frequency detector circuit (5) comprising a 
time-axis to frequency-axis converter circuit (2) 
for converting an input data stream that is 
obtained by sampling a received signal on a 
time axis into a power spectrum on a frequency 
axis, and a MAX search circuit (3) for search- 
ing, in the observation frequency range deter- 
mined by a sampling frequency, for the 
frequency corresponding to the maximum 
value of the power spectrum, 
a frequency computing circuit (60) for estimat- 
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